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Abstract. The flow generated by low-density radiatively 
driven winds which decouple their gas and radiation fields 
is discussed. In particular we concentrate on flow which 
is still bound to the star and can therefore reaccrete. The 
wind decelerates after decoupling and eventually stalls. A 
shell of gas is generated, and we find that this shell is 
unstable and contracts back to the star with periods of 
hours to days. 

We find that the pulsating shells may be difficult to 
observe, as their emission is variable and the maximum 
emission at H Q (of ~ 1% of the continuum) occurs over a 
small fraction of the shell cycle. 

Key words: stars: mass-loss - hydrodynamics - stars: 
early type - circumstellar matter 



1. Introduction 

The winds from hot stars are thought to be generated by 
radiation pressure on optically thick UV resonance lines, 
and the theory of line-driven flow is very successful in 
accounting many observed wind features (see Castor, Ab- 
bott & Klein 1975, Abbott 1980, Pauldrach et al. 1986, 
Kudritzki et al. 1989). 

In low density winds, however, there exists the pos- 
sibility that the radiation force and the wind flow may 
decouple (Springmann & Pauldrach 1992, Porter & Drew 
1995 [hereafter PD95], Babel 1996). The decoupling pro- 
cess strips the metallic ions from the rest of the plasma, 
and as the radiative force on the flow is mediated by the 
ions, then the wind receives no further acceleration. The 
winds which are most likely to undergo this decoupling 
are B star winds (Babel 1996, PD95), and metallic A star 
winds (Babel 1995). The frictional interaction between the 
metallic ions and the rest of the wind may also seriously 
interfere with the radiative equilibrium (Springmann & 
Pauldrach 1992, Gayley & Owocki 1995). Decoupling radii 
for low density winds may be close to the photosphere for 
B stars (PD95) or indeed may be associated with the pho- 
tosphere on the case of some A star winds (Babel 1995), 
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where there is no region outside the photosphere where a 
fully coupled wind exists. 

The line-driven wind accelerates normally when fully 
coupled to the radiation field, but once decoupled cannot 
receive any further acceleration. It is possible that the 
wind will decouple before it has reached escape velocity, 
in which case the decoupled flow will stall at some radius 
and fall back toward the star. It is this aspect of decoupled 
flows which is examined in this paper. 

The rest of this paper is structured as follows: in §2 the 
physics of wind decoupling is examined, and in §3 hydro- 
dynamical simulations of decoupled winds are presented. 
The observational signatures of the shells are presented in 
§4 and a discussion and conclusions are given in §5. 

2. Decoupling physics 

There are two ways in which the radiation and matter 
fields in a line-driven wind may decouple: ion stripping 
(PD95, Springmann & Pauldrach 1992) and shock decou- 
pling (PD95, Krolik & Raymond 1985). 

The physics of ion stripping was first noted in the con- 
text of electrical conductivity by Dreicer (1959, 1960) and 
has its roots in the basic mechanism allowing optically 
thick metallic ion lines to mediate the force on a wind. 
The metallic ions are accelerated via photon scattering off 
their UV lines. The ions then share this acceleration with 
the rest of the wind (hydrogen and helium ions) by a pro- 
cess similar to friction. This frictional interaction depends 
on the relative drift velocity of the ions through the rest 
of the wind. For low drift velocities the force is propor- 
tional to the drift velocity, however it reaches a maximum 
when the drift and thermal kinetic energies are equal. Be- 
yond this the frictional interaction decreases rapidly with 
increasing drift velocity. Therefore if the ions' drift veloc- 
ity does become large enough, then as there is little or no 
frictional interaction with the rest of the wind, they may 
freely accelerate. This leaves the rest of the wind with no 
acceleration, and it will then just be acted upon by the 
gravitational attraction of the central star. 

Shock decoupling occurs in a different way. It has been 
known that hot star winds may be unstable to the growth 
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of instabilities for some time (Lucy & Solomon 1970, Mac- 
Gregor et al. 1979, Owocki & Rybicki 1984) although it is 
unclear whether the winds are inherently unstable or only 
advectively unstable (see the review by Owocki 1991). If a 
low density wind passes through a strong shock, then the 
postshock gas may be too diffuse to cool radiatively, and 
so the gas remains in its high temperature, superionized 
state (Krolik & Raymond 1985). The wind may decou- 
ple from the radiation if the postshock gas remains out of 
equilibrium, as the main metallic species responsible for 
line-driving have been ionized away (Castor 1987). 

Although both of these processes may in principle oc- 
cur, only the former (ion stripping) is considered here as 
it is more deterministic - the decoupling radius for shock 
decoupling can only be calculated if the star's X-ray emis- 
sion (to determine the shock velocity) and mass-loss rate 
are available, whereas for ion stripping, only the mass-loss 
rate is required. It should be noted that the decoupling ra- 
dius is similar in each case (see PD95). 

For ion stripping, the decoupling radius Td may be es- 
timated by considering the relative velocity between the 
ions and the rest of the wind. As it is this drift velocity 
which sets the frictional interaction then it also defines 
how much acceleration is imparted to the rest of the wind. 
When the drift velocity becomes comparable to the ions 
thermal speed, the frictional interaction reaches a maxi- 
mum beyond which it falls with increasing drift velocity. 
At this point the frictional interaction between the ions 
and the rest of the wind becomes small and the ions are 
(nearly) free to accelerate out of the wind. This radius is 
where the matter and radiation fields effectively decouple. 
For a beta-velocity law this yields 

M_ 9 1 

(see Springmann & Pauldrach 1992 and also Gayley & 
Owocki 1995) where Voo 3 is the terminal velocity of the 
wind measured in units of 10 3 km s _1 , and the stellar ra- 
dius i?* is expressed in solar units. Also Yi is the mass 
fraction of the ions, Zi is the degree of ionization, and 
M-g = M/10 _9 M Q yr _1 is the mass loss rate derived 
from UV line profiles. T4 is the temperature of the gas 
in units of 10 4 K and InA is the Coulomb logarithm. The 
approach of calculating the ion-stripping radius used to 
derive eq.l is, however, strictly only valid when the wind 
is fully coupled. 

Therefore, a more exact calculation needs to be under- 
taken including the precise frictional force expression (e.g. 
see Springmann & Pauldrach 1992, Dreicer 1959, 1960). 
Although the frictional heating may increase the temper- 
ature of the wind and cause decoupling at smaller radii 
(PD95 noted this fact), the decoupling radii and velocities 
for a model B2 star have been calculated assuming isother- 
mality. Temperature profiles of coupled radiatively-driven 
winds have been calculated by Drew (1989) using a full 
description of radiative heating and cooling processes in 
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Fig. 1. Decoupling radii (left) and velocity (right) for a 
7.5M Q , A.ORq , 2xl0 4 K star as a function of mass-loss 
rate. The escape velocity is 840 km s -1 and the wind ter- 
minal velocity is 2260 km s _1 . 

the wind. Drew's calculations show that the wind temper- 
ature falls below 0.6T o ff only for r > 2i?» , although the 
frictional heating from the relative drift of ions through 
the rest of the wind was not included. Therefore the as- 
sumptions made here that the wind temperature is 0.8T c ff 
(Klein & Castor 1978), and that the flow is isothermal is 
probably not seriously in error. 

For all the calculations undertaken in this paper, we 
use a "standard" B2V star which has stellar parameters of 
M* = 7.5M , R* = A.0R Q , T cfr = 20,000K. The escape 
velocity of this star is v csc = 840 km s" 1 , and the fully 
coupled terminal wind velocity is Voo = 2.2av csc /(l — a) = 
2260 km s" 1 (Friend & Abbott 1986), for the line-driving 
parameter a = 0.55 (see Castor, Abbott & Klein 1975). 
The decoupling radii and velocities have been calculated 
for our standard star, and are presented as a function of 
mass-loss rate in fig.l. 

The escape velocity of the star is 840 km s" 1 and so 
from fig.l, it can be seen that if the mass-loss rate is less 
than - 10- 9 - 2 M Q yr" 1 (= 6.3 x 10- 10 M Q yr" 1 ), then the 
wind will decouple and still be bound to the star. 

3. Results 

3.1. Analytic considerations 

If the radiation and matter fields decouple at a velocity 
Vd < v csc then, as noted above, the material is still bound 
to the star. In this case the local radial velocity must be 
calculated for flow in a gravitational field. The radial ve- 
locity will decrease until the flow stalls at an outer radius 
r s which can be estimated by equating the change in po- 
tential and kinetic energies, leading to 

r s _ [R, 1 /i?,fl 2 \ i^A 2 !" 1 (2] 
R* ~ [r d 2 \GM* J \ a) \ ' { ' 

As there is flow still with positive radial momentum be- 
hind this stalling front, a shell of stalled wind will form 
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being supported by the ram pressure of the gas behind it. 
If we assume that the shell can be supported by the wind 
at some radius r s h e ii, then we can estimate the maximum 
total mass in the shell. The supporting ram-pressure P ram 
of the stalling gas will be pv 2 where p is the density and v is 
the local radial velocity The oppositely directed pressure 
generated by the shell is the due to the attractive grav- 
itational force of the star GM* M s h e ii/''fheii where M srie ii 
is the shell mass and G is the gravitational constant. This 
force acts over the surface area of the shell 47rr 2 hell . Bal- 
ancing the two oppositely directed pressures yields 



M 



Mvr; 



shell 



hell 



GM, 



(3) 



Here we have used the mass-continuity equation for spher- 
ical flows M = 4irr 2 pv, where M is the mass-loss 
rate of the wind. If v 2 = f/100 km s _1 and M-g = 
M/lO^Moyr- 1 , then 



^ « 1.2 x 10- 14 M_ 9U2 ( Ji- 



Mo 



(4) 



Note that a distinction has been made between the stalling 
radius r s and the radius of the shell r s h c ii- This is because 
the shell will only appear at the stalling radius if the wind 
first expands out into a vacuum. If there is any ambient 
density around the star p am b, then this is swept up into 
the shell when the wind is "turned on" . 

The evolution of the wind will broadly be composed 
of several parts. Initially the wind will be driven into the 
surrounding medium and will sweep up a shell of ambient 
gas. This shell will stop expanding when the shell and ram 
pressures balance. The shell then grows as it is receiving 
gas from the wind (at presumably the same rate as the 
mass- loss rate from the star). Finally the shell becomes 
too massive to be supported and will fall back toward the 
star. 

If we assume that p am b is small then the shell ra- 
dius is the stalling radius. Also the wind's ram pressure 
in this case is the thermal pressure pa 2 , where a is the 
sound speed. For stalling radii of r s ss 5i?* = 2Oi? 
the shell mass for our standard star with a mass-loss rate 
of 10 _9 5 Mq yr _1 which will be supported by the stalling 
wind is M s h c ii ~3x 1CP 14 M0 . Given that the shell mass 
grows at a rate of M, we find that it is stable for a time 



t = 380v 2 



R, 



Mr, 



(5) 



which for the example above is t « lhr. The timescales 
here are similar to the flow timescale of the wind (approx- 
imately i?„ jv w lhr for a wind speed of 1000 km s _1 ). 

After this time the shell's mass will have increased such 
that there is insufficient ram pressure to levitate the shell 
and so it will fall back to the star. As the shell falls inward, 
the wind ram pressure increases - both the wind density 



and the radial velocity v increase (Mach numbers in radia- 
tively driven flows can be ~ 100). However, the increasing 
mass of the shell is now acting over a smaller area. Hence 
the effective pressure exerted by the shell inward due to 
the gravitational attraction by the central star increases. 
This pressure of the shell is the gravitational force divided 
by the area of the shell: P s hcii oc r~ 4 whereas the ram pres- 
sure is Pram oc v/r 2 . The ratio of these two pressures is 

-2 



Phdi 8.6 x 10" 



AL 



shell 



, . • (6) 

M^v 2 \M Q J \ M Q J \R Q ) 

Note that in the decoupled part of the wind v <~ r -1 / 2 , 
and so as M s h c ii increases and moves toward the star (r 
decreases) the ratio above grows. This then implies that 
the shell can not be stable (balanced) at any radius, and 
must fall back to the star. 

Once the shell has contracted to the star, the wind is 
free again to expand to its stalling point and the cycle 
restarts. This analysis predicts that stars with bound de- 
coupled winds will inevitably produce periodic structures 
associated with the wind. However, a word of caution is 
warranted here : this scenario predicts that a shell of gas 
should form and collapse onto the star with timescales of 
hours-days. There is a possibility that the oppositely di- 
rected gravity and pressure gradients as the shell contracts 
may break up the shell due to Rayleigh- Taylor instabili- 
ties. 

3.2. Numerical hydrodynamic modelling 

To illustrate the preceeding section's scenario, a ID nu- 
merical hydrodynamic simulation has been completed. 
The computer code utilises the second order Godonuv 
scheme due to Falle (1991) to solve the hydrodynamic 
equations. The line force is calculated using the stan- 
dard Castor, Abbott & Klein (1976) formalism of the 
force multiplier, supplemented with the finite-disc cor- 
rection (Friend & Abbott 1986) and the density weight- 
ing factor from Abbott (1982). The grid has a spacing 
of 0.005-R* at r =R* , increasing by 1.5% to a maximum 
spacing of 0.09P* at the outer radius of 6.75P* (there are 
200 grid points). The outer boundary allows gas to flow 
through it freely, whereas the inner boundary has the den- 
sity fixed at 10~ 12 g cm -3 , and the velocity extrapolated 
from the neighbouring zones. We can expect that after a 
time the radial velocity at the outer boundary will become 
negative and matter will flow into the computational do- 
main. This causes the density to rise at the outer parts 
of the simulation, and finally lead to a significant ram 
pressure acting on the shell from infalling material. Ulti- 
mately this leads to the central star being surrounded by 
high density, infalling material, completely thwarting even 
the initiation of the wind! Clearly this is an artifact of the 
outer boundary. After experimentation, it is found that 
the outer boundary which illustrates the scenario most ef- 
fectively is a "valve", i.e. the boundary is free-flow when 
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the radial velocity is greater than zero, but otherwise we 
set the radial velocity to zero, and the density to a small 
value (10~ 30 g cm -3 ). 

In order to examine the effect of differing mass-loss 
rates on the flow we prepared five fully-coupled wind so- 
lutions as initial conditions. These have been generated 
for our standard B2 star. We find that the mass-loss rate 
for this star is ~ 10 _8 ' 3 Mq yr _1 . As we require differing 
(and much lower) mass-loss rates to assess its effect on 
the decoupled flows, we have reduced the radiation force 
by as much as an order of magnitude. This may be seen 
as a change in metallicity - the radiative acceleration is 
proportional to the metallicity raised to the power of 1 — a 
(Abbott 1982, Castor, Abbott & Klein 1975). The lowest 
metallicity model implied is then 0.006 solar. By reducing 
the line force we have generated simulations with have 
mass-loss rates in the range (1.6-4.4) x 1O _1O M0 yr _1 (see 
table 1). We have fitted a "beta" velocity law to the first 50 
points (upto 1.3i?* ) of the initial conditions and find that 
they are all well fit by v =2100(1 - /r) - 8 km s"\ 
although for the models with higher mass loss rate, /3 falls 
to 0.73. With this velocity structure, we calculate the de- 
coupling velocities and radii as described in §2. These are 
presented in table 1. Note, we have assumed that all mod- 
els have solar metallicity for the calculation of the decou- 
pling radii, which may be inconsistent with the generation 
of the initial conditions if the low radiative acceleration is 
interpreted as due to a change in metallicity. 

In order to simulate the decoupled flow, we use these 
fully coupled wind structures as initial conditions and sim- 
ply turn off the radiation force at the decoupling radius. 
The flow is calculated in each simulation for at least ten 
shell episodes to ensure that the initial conditions have no 
effect on the subsequent flow. 

Part of the simulation's results for model A are shown 
in fig. 2 - the panels of velocity (left) and density (right) 
are labelled with the time in hours (the starting point for 
the time has been shifted to coincide with fig. 3). All the 
models have qualitatively the same behaviour with only 
the physical scales changing from model to model. The 
first couple of cycles have not been shown to remove any 
start-up transients. The dotted lines on the panels are the 
starting condition described above as a reference. The first 
panel shows the dense shell forming at a radius of 3.5i?* . 
The shell's mass increases (to ~ 1.8 x 10~ 13 M Q ) until its 
weight is large enough to overcome the wind ram-pressure, 
and so it contracts toward the star (at around 5hrs later 
- the second panel of fig. 2). The shell gains mass and falls 
faster with time (panels 3 and 4), finally collapsing onto 
the star just before 37hrs. Once the shell makes contact 
with the star, the wind is again free to accelerate to de- 
coupling. This process sweeps up some of the material 
left behind by the previous cycle's contracting shell and 
a new shell is generated moving outwards (the lower two 
panels of fig. 2). The whole cycle takes around 28hrs. Al- 
though the timescale for the shell growth is very similar 




radius (R + ) radius (R^) 

Fig. 2. Numerical hydrodynamic modelling of a the for- 
mation and evolution of a spherical shell (Model E). The 
dotted line is the steady, fully-coupled starting condition. 
The plots are labelled in time (in hours). The radiation 
force was set to zero at radii larger than 1.27.R* . 

to our previous estimate, the timescale for the whole cycle 
is ~ lday for this star. The periods for the shell episodes 
for the different models are shown in column 7 of table 1. 

We have calculated the emission measure EM= 
J n 2 e dV (where n e is the electron density, and V is the 
volume of the shell) for the simulations. The median max- 
imum value for the emission measure from all the simu- 
lated cycles is shown in column 6. Fig. 3 shows how the 
radius, velocity, and emission measure changes with time 
through a shell episode. We have identified the shell at the 
point where there is a local maximum in density (when 
the shell is close to the star it becomes difficult to identify 
unambiguously and so these points have been left out of 
fig. 3). The wind initially sweeps up the ambient gas cre- 
ating a low density-contrast shell. In this phase the shell 
is decelerating from the decoupling radius, and moving 
out in radius. Finally the shell comes to rest (v s = at 
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Table 1. Model parameters for the simulations, and Vd are the decoupling radius and velocity respectively. The 
stalling radius r s is calculated from eq.2. The shell period and the median maximum emission measure (over several 
cycles) are listed in columns 6 and 7. 



Model logioM (M Q yr 1 ) r d (R* ) v d ( km s 1 ) r s (R« ) logipEM period (hr) 



A 


-9.35 


1.27 


590 


3.50 


56.0 


27.8 


B 


-9.46 


1.26 


568 


3.09 


55.7 


20.8 


C 


-9.59 


1.23 


522 


2.39 


55.1 


14.6 


D 


-9.68 


1.20 


447 


1.82 


54.6 


8.9 


E 


-9.78 


1.17 


378 


1.53 


54.1 


6.6 



71 



r w 3.5i?» ), and it starts to gain mass. The emission mea- 
sure starts to increase significantly during this time. When 
it cannot be supported any longer the shell starts to fall 
back toward the star, finally coming into contact with the 
star at a velocity of around -500 km s _1 . 

All of the simulations produce very similar plots except 
the scaling in time and position are different. We find that 
the stalling radii are very similar to those predicted using 
eq.2 (table 1, column 3). 

The period for initial shell growth has already been 
shown to be dependent on the mass-loss rate via the 
stalling radius (eq.5) which is, in turn, set by the de- 
coupling radius. As the mass-loss rate increases, the de- 
coupling velocity/radius increases (see table 1), and so 
the stagnation radius (at which the shell forms) also in- £ 
creases. We expect that the scaling relationships in the ^ 
shell-growth timescale (eq.5) is similar for the timescale 
of the shell episode as a whole, although the time for shell 
growth is only a fraction of the shell cycle. Hence we ex- 
pect that the period should be proportional to the stalling 
radius squared. In fig. 4 the logarithmic period-stalling ra- 
dius relationship is displayed. We have performed a least- 
squares fit to the above data and find that the period P 
is related to the stalling radius via P = 3.22rJ- 7 which is 
close to the expected scaling. We also find that the max- 
imum emission measure follows the same scaling with ra- 
dius. 



4. Observable signatures 

The total mass in the shell will typically be the wind 
mass-loss rate multiplied by the total lifetime of the shell. 
Timescales of ~day have been found from the previous 
section, yielding maximum shell masses of M s hoii ~ 10~ 13 
10 _12 M© . Clearly such a low mass will not have a large 
influence on the gross observational properties of the star. 
The most obvious effect of decoupling is that the max- 
imum wind velocity will be observed to be lower than 
the terminal velocities from theoretical studies - a phe- 
nomenon which has been known observationally for some 
time (Grady et al. 1987). 

Would the shell give rise to any variation in the lines 
of the star? Clearly there would be no change in the UV 
line profiles until the shell radius becomes less than the 
decoupling radius. Here the maximum edge velocity will 




20 40 
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60 



Fig. 3. Time resolved attributes of the shell for Model A. 
Position (top) and velocity (middle) of the shell in time 
through a shell episode. The lower panel shows the total 
emission measure of the wind structure. The solid dots 
indicate the times of the corresponding panels in fig. 2. 



decrease as the shell falls toward the star. However, the 
occurs only over a small fraction of the time of the whole 
cycle. 

It may be possible that the shell may provide a vary- 
ing component to the hydrogen lines. Let us assume that 
the shell is always optically thin and calculate the ex- 
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Fig. 4. The log-log relationship between the period of shell 
episodes P and the stalling radius of the wind r s from 
eq.2. The straight line is the least-squares fit to the points: 
P = 3.22(r s /R* ) 17 hours. 

cess emission at H Q . This is simply j = hvoF§ji e n p j 'Air 
erg cm~ 3 s _1 st _1 where n e and n p are the electron and 
proton number densities, and a c ^ a is the Ha recombina- 
tion coefficient of 5.96 x 10~ 14 cm 3 s _1 assuming Case B 
recombination at 20,000K (Osterbrock, f989). Assuming 
that the emission is optically thin then the shell produces 
a flat-topped emission profile with a range in wavelength 
given by the velocity of the shell. When the shell is close 
to the star, then the star occults the far side of the shell 
and so curtails the maximum blucshiftcd emission. If the 
angle between the pole and a point on the shell is 6, then 
emission from that point is Doppler shifted by 



SX = X — sin 6* 

c 



(7) 



where v s is the velocity of the shell, and c is the speed 
of light. Combining this with the expression for j, and 
integrating the volume emitting at the same wavelength, 
we derive the emission as 



hcttf n e n p 



r 2 dr ( — I erg s x st 1 hz 



(8) 



We assume that that star is emitting as a black-body 
and calculate the fraction of the black body flux the shell 
emits (of course in a real star there is also photospheric 
absorption) . After some simple manipulation we find that 
this ratio is 



K = 



hv/kT 



1) 



c° a 



EM 



is 3 AirR* 2 v s 



(9) 
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Fig. 5. Grey-scale representation of the time series of the 
excess line emission at H a due to the shell for Model A 
(see text). The grey scale in linear, and the darkest tones 
correspond to f % of the continuum emission. 



where the symbols have their usual meanings. This extra 
component extends over a wavelength range 



Ao 1- 



r s 



< X < 



(10) 



where Ao = 6563A is the central wavelength of the line. 

We have calculated the extra component of H Q for 
Model A and present it as a time series in fig. 5. The grey 
scale is linear with a range of 0-1% of the continuum with 
the the darkest tone corresponding to 1%. We find that 
the maximum value of 1Z for Model A is ~ 1%, i.e. the 
flat-topped shell component has a contrast of ~ 1% of 
the continuum flux. Whilst this is clearly an observable 
amount, it is only present at significant values for a small 
of the cycle. 
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5. Discussion & Conclusions 

We have demonstrated that low density radiation-driven 
flows may decouple before the wind becomes unbound 
from the star, and therefore will reaccrete. Abbott & 
Friend (1989) discussed a model of a line-driven wind with 
a line force cut-off to mimic the effects of shocks in the 
wind. However, their models had force cut-offs such that 
the wind beyond decoupling was unbound and therefore 
would escape from the star. 

This is the first suggestion that these flows may still 
be bound to the star post-decoupling. We have examined 
the attributes of the periodic shell structures formed in 
the flow and have calculated the line emission. 

From our analysis we conclude that these shells may be 
observable in a line monitoring campaign, with the largest 
emission coming from stars with winds which decouple 
close to the escape velocity. These stars will also have the 
longest period for the shells making it more probable that 
the shell can be observed. 

It is worthy of note that this process of reaccretion will 
change the metallicity in the outer layers of the star over 
a long time period. As it is only the metals which become 
unbound, then their abundance decreases with time in the 
outer layers when the shells reaccrete. In a future study 
we will address the case of rotating winds which may be 
applicable to classical Be stars as these stars do seem to 
have small mass-loss rates. 
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